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In this article we show that the inhomogeneous density obtained from a density-functional theory 
of classical fluids in the canonical ensemble (CE), recently presented by White et al [Phys. Rev. 
Lett. 84, 1220 (2000)], is equivalent to first order to the result of the series expansion of the CE 
inhomogeneous density introduced by Gonzalez et al [Phys. Rev. Lett. 79, 2466 (1997)]. 
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A statistical mechanics ensemble is a collection of identical systems under the same external conditions. Although 
the choice of a particular ensemble for studying a concrete system should be guided by the conditions in which the 
system is found, one can choose — due to mathematical or computational convenience — any ensemble for analyzing 
the equilibrium properties of the system. This way of proceeding, based on the equivalence of the ensembles in the 
thermodynamic limit, is only justified for systems with a very large number of particles. For small systems, however, 
the ensembles are no longer equivalent and the external conditions must determine the choice of ensemble. 

In this context, the use of density-functional theory (DFT) for the study of classical inhomogeneous fluids has 
been usually limited to the grand canonical ensemble (GCE), where the temperature T and the chemical potential 
yU are fixed by an external reservoir. A large variety of inhomogeneous situations has been successfully studied by 
means of DFT in the GCE |l|-^ . These situations include fluids confined in narrow pores or capillaries |Q , or even 
spherical cavities P-Rj, which are implicitly assumed to be open, i.e., allowing exchange of particles with a reservoir. 
This assumption is crucial for situations with a small number of particles where, depending on the choice of ensemble, 
important differences may arise in the equilibrium microscopic structure of the system [nj|] . If one wishes to investigate 
the properties of a small closed system at temperature T, the study must be performed in such a way that one obtains 
results in the canonical ensemble (CE) because the number of particles N is fixed. In DFT this goal can be achieved 
by means of two different approaches. On one hand the DFT could be formulated in the canonical ensemble ^, with 
a minimum free-energy principle with fixed T and N , and an appropriate CE functional. Very recently, this approach 
has been explicitly realized lln] by considering an approximate expression for the CE functional. On the other hand, 
one can perform a conventional DFT study in the GCE and then relate the obtained properties to those of the CE. 
This approach was followed in Refs. |R,|[ where the CE density profile of a hard-sphere fluid in a small spherical cavity 
was calculated by means of a series expansion in terms of the corresponding GCE profile. The aim of the present 
paper is to show that these two approaches yield equivalent results to order \/N . For clarity we start with a brief 
summary of the main results of the two approaches. 

The first approach is based on the following series expansion of the CE density profile pc (r) in terms of its corre- 
sponding GCE density profile Pgc(r): 

Pe(r) = pgc(r) - \^\N)^p,^{v) + O(^) (1) 

where the grand canonical profile is obtained for a chemical potential fi that leads to an average number of particles (N) 
equal to N — the fixed integer number of particles in the CE, and A^(A^) = {N'^) — (N)'^ is the mean square fluctuation 
of the number of particles in the GCE. Higher-order terms in the above expansion also depend on fluctuations and 
on variations of the GCE profile w.r.t. (N) mH- In DFT, the grand canonical profile pgc is the solution of the usual 
GCE Euler-Lagrange equation with chemical potential /i and external potential 14xt(r), 
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where .^gc[p] is the GCE free-energy functional. For a chemical potential leading to a given {N), Eq. (||) can be 
rewritten as M: 



Pgc(r) = (A^)exp [-/JFextW +c(i)(r; [pg,])] / J drexp [-f3V,^t{r) + c^^\r; [pge])J , (3) 

where f3 — 1/k^T is the inverse temperature and c-^> is the one-body direct correlation function 



cW(r;[p]) = -/3 
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being J^gc-id the usual ideal-gas free-energy. This correlation function is the first member of the direct correlation 
hierarchy 



c(")(ri,...,r„;[p]) = -/3 
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From Eq. (Q), in DFT it is possible to obtain density profiles normalized for a given {N). This allows for obtaining 
approximate CE profiles using (m where the derivatives w.r.t. {N) are calculated numerically. (We note that, using 
the thermodynamic identity A^(iV) — d{N)/d{/3iJ,), the mean square fluctuation can also be expressed as a derivative 
w.r.t. {N}.) This procedure was used in mS to obtain the CE density profile of a hard-sphere fluid confined in a 
hard spherical cavity. 

The second approach consists of an approximate expression for the free energy functional in the CE. On the basis 
of the standard saddle point relation between the CE Helmholtz free energy and the GCE grand potential [Q the 
following approximation for the CE free-energy functional jFc was proposed in Ref. llfl] : 



PTM « /3^ge[p] + 2 log 2^ A2 (TV; [p]) 
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where the functional dependence of the GCE mean square fluctuation A^(A^) is made explicit. Since we are now 
working in the canonical ensemble, the equilibrium density profile Pc(r) is obtained by minimizing the functional 
!Fc[p] + J drp(r)V;,xt(r) subject to the constraint 



/ drpc(r) = N. 
Using the Lagrange multiplier technique one obtains pl]| 
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+ Kxt(r) = A 
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where the Lagrange multiplier A must be calculated from the constraint (M). This equation can be re-expressed as 



Pe(r) = A^exp[-/3Kxt(r) +c(i)(r; [pc]) + ^(r; [pc]) 



dr 
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-/3yext(r) + c(i)(r;[pe])+e(r;k]) 



where 
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which, for the saddle-point (SP) approximation (|6|), becomes 
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In order to calculate Csp(r; [p]) it is important to express the mean square fluctuation A^(A^; [p]) as a functional of 
the density. This can be done conveniently in the GCE by means of the density-density correlation function 0,|l^ 



a(ri,r2)=r'- 
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^(A*-Kxt(r2))' 
since this function normalizes to the mean square fluctuation, i.e., 



A\N) ^ j j dv^dv2Q{vi,V2) . 
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In addition, taking into account that Q is the functional inverse of the second derivative of the GCE free-energy 
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^ <5(ri-r2)-c(2)(ri,r2), (14) 
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and satisfies the Ornstein-Zernike relation P,p^ 

f dr2g-\ri,r2)g{r2,r3) = 5(ri - rg) , (15) 



one obtains |1C 



A\N;[p]) = Jdrr{r) (16) 

^ -A^N;[p]) = [ /drirfr2^^-^i^r(ri)r(r2) (17) 
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/ / dridr2c(3) (r, n, r2)r(ri)r(r2) (18) 



where 



r{r) = J dr,g{r,r,) , (19) 

is obtained from the following averaged Ornstein-Zernike relation 

r(r) = p{r) + p{r) J drir(ri)c(=^)(r, n) . (20) 

In deriving (n^ we have considered the functional derivative w.r.t. density of the Ornstein-Zernike relation (nsh and 
exploited the fact that Q and Q~^ are functional inverses. We note that the key difference between the GCE result 
(^ and the CE density (^) is the term ^(r; [pc])- We also note that, in this approach, using (||) one directly obtains 
the CE profile while, in the previous approach, the result (||) of GCE-DFT had to be inserted into (y) in order to 
obtain an approximation for the CE equilibrium density. In what follows we shall show that both approaches, agree 
to first order, though they yield different results due to higher-order contributions in the saddle point approach. We 
first derive some useful relations and then we show the equivalence to first order of the approaches. 

Our starting point is the well-known result of GCE density functional theory that, for given intermolecular potential, 
chemical potential /i, and temperature T, only one external potential can determine a specified equilibrium density 
profile jljQi. Thus there must exist an external potential t4xt(r) so that its corresponding GCE equilibrium density 
/9gc(r; [Kixt]) (here the functional dependence of pgc is made explicit) is equal to the CE result pc- Performing a 
functional expansion of /Ogc(r; [Kjxt]) about T4xt and using definition ( p^ ) we obtain: 



/?-iAp(r) =/'dri5(r,ri)AT/cxt(ri) 



where 



and 
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Ap(r) = pgc(r; [Kxt]) - Pgc(r; [Vcxt]) = Pc(r) - Pgc(r) (22) 



A14xt (r) = Fcxt (r) - Fext (r) . (23) 

Therefore, Eq. (21) provides a link between Ap and AV^xt via a functional expansion where the coefficients belong 



to the standard distribution function hierarchy. 



At this point we would like to emphasize the role played by t4xt in the present work. t4xt is the external potential 
that, at chemical potential /x, yields the canonical profile pc in a GCE approach and this implies that p^ is the solution 
of the following GCE Euler-Lagrange equation 



(5J^gc[p] 
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+ Kxt(r)=/i- (24) 



This fact makes meaningful the use of functionals of pc like Q ^(ri,r2; [pc]) — f3S{p — Vc:^t{^i))/Spc(j^2), its inverse, 
g{ri,r2;[pc]) = G{ri,r2;[Vcxt[pc]]) = P^^dpc{ri)/S{p - Vcxt{iC2)), or the mean square fluctuation A'^ {N ; [pc]) . On 
the other hand, since /Oc(r) is the equilibrium density profile in the canonical ensemble, it is the solution of the CE 
Euler-Lagrange equation (ra) where, in comparison with this GCE equation, the free energy is Tc — the CE functional, 
the external potential is V^xt, and the Lagrange multiplier A is used in the place of the chemical potential p. From 
Eqs. (^) and (|2j) and definition (|l^) we obtain 

Ayext(r) = A-/i + /3-iC(r;[Pc]). (25) 

In the uniform limit, where Pgc(r) — > Po = {^) /^ — ^1^ and also Pc(r) ~^ Po, from (||) and ( p^ one has AVexti'r) = 
and thus 

p-X^p-'aPo), (26) 

and ( |25| ) can be rewritten as 

/3AKxt(r)=e(r;[Pc])-e(Po). (27) 



Using the saddle-point approximation ^gp [Eq. (pi])], this expression could be employed in (21) to obtain an approx- 
imation for the difference Ap. Conversely, since Eq. (^ gives an approximation for Ap, an expansion inverse to 
( pl| ) would provide a way to obtain ^. This inverse expression can be easily derived by substituting the functional 
expansion of SJ^gc[p]/Sp{r) about pgc in Eq. (P4|). Using definition ([l4|), we obtain 



/?AUext(r)= f drig-\r,ri)Ap{ri) 



+ 1 /dridr2^^-^^Ap(ri)Ap(r2) + ... (28) 
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where we have exploited the fact that Pgc(r) is the solution of the usual GCE Euler-Lagrange equation (g). 

Expansions ( pi| ) and (Eq) are asymptotically exact relations linking Ap and AV^^t- However, these expansions need 
to be truncated in order to become suitable for practical applications. In particular, to first order, Eq. ( pi|) becomes 

/3-iAp(r) « /"drig(r,ri)AUext(ri), (29) 

and Eq. (28) reduces to 

/3AUext(r) «ydria-i(r,ri)Ap(ri). (30) 

Approximations ( p9|) and ( pOJ ) are, by virtue of Eq. (|l5|), equivalent equations; this fact shows the consistency of the 
truncation of the expansions. Either ( p9| ) or ( pOJ ) provide a simple (first order) relation between the differences Ap 
and AT4xt- Using these equations, we shall show that the approximation (|l|) for Ap is equivalent to first order to the 
saddle point approximation ^sp [Eq. (pi])]. By considering the derivative of the GCE Euler-Lagrange equation (0) 
w.r.t. (N) at fixed Kixt, one obtains the exact relation 
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where the chain rule for functional differentiation together with definition (|lj) and the identity d{N)/d{(3p) — A'^{N) 
have been used. This equation can be rewritten, via Eqs. (p5[) and (p9|), as 



d{N) A^{N)' 
Differentiating this equation w.r.t. (N) and using Eq. (yj), we obtain 
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which inserted into ( |30| ) yields the following approximation for Al^xt 

.2 



where we have used the identity 

/"drig-i(r,ri)r(ri) = l (35) 

which follows from Eqs. ( JSl] ) and (^) [or, cquivalently, from ( p^ and (|9|)]. Considering the derivative of Eq. ( |35| ) 
w.r.t. (A^), and using the chain rule and (P2|), Eq. (gj) can be re-expressed as 

1 fdA^N) SA^N)\ ^3^^ 



2A2(7V) \^ d{N) <5pgc(r) ^ ' 

where in the last equality we have used (^3)- Comparing (p7| ) with (g^ we obtain 



1 5A^{N) 
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and 
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where ^gp was defined in (|ll|). Finally, we note that in (38) it is shown that S,[pc\ is approximately equal to ^sp[Pgc], 



i.e., ^sp evaluated at pgc instead of pc- This approximate equality also holds for ^sp[pc] as we show in what follows. 
Expanding ^gp about pgc we obtain 

U(r; [pc]) = e.p(r; [Pgc]) + / dri^^^lMlAp + • ■ • (40) 

and, taking into account that ^ is already a 0(Ap) quantity, 

e.p(r; [pe]) « esp(r; [pgc]) + 0(Ap)2 (41) 

This proves the equivalence to first order of the two approaches for obtaining an approximate density profile in the 
canonical ensemble. 

In summary, we have shown that two different approaches for obtaining the density profile of a fluid in the canonical 
ensemble are equivalent to first order. The demonstration was based on considering an external potential V^xt for 
which the equilibrium density in the grand canonical ensemble is precisely the canonical ensemble result. Using this 
external potential we have been able to work in the framework of the grand canonical ensemble where approximations 
similar to those carried out in this work are commonly encountered. 

The proof of the equivalence gives additional support to the saddle point approximation for the CE free-energy 
functional introduced in |jlffl . This approximation allows for a CE-DFT treatment of fluids conflned in a closed cavity 
with excellent agreement with simulation data. However, the SP free-energy functional was proposed on the basis of 
the well-known SP relation between the equilibrium CE free-energy and the grand potential of a homogeneous fluid, 
assuming that this relation would be also a good approximation for inhomogeneous fluids. This assumption is thus 



reinforced by the results of the present paper which are vahd for any inhomogeneous situation. As a final remark, we 
would like to stress the fact that our demonstration has focused on the approximate CE density rather than on the 
free-energy functionals, and the equivalence between the two approaches must be understood in this sense. 
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